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Abstract. Alfven waves are considered as potential transporters of energy 

^^ ' heating the solar corona. We seek spectroscopic signatures of the Alfven waves 

^^ . in the chromosphere occupied by a bright network element, investigating tem- 

• ' poral variations of the spectral width, intensity, Dopplershift, and the asym- 

i-i^ I metry of the core of the Ha spectral line observed by the tunable Lyot filter 

S-^Hi installed on the Dutch Open Telescope. The spectral characteristics are de- 

Q ' rived through the fitting of five intensity samples, separated from each other 

$-H , by 0.35 A, by a 4*'^-order polynomial. The bright network element displays the 

[/3 ■ most pronounced variations of the Dopplershift varying from to 4kms~^ 

C^ ' about the average of 1.5kms~^. This fact implies a persistent redshift of the 

Ha core with a redward asymmetry of about 0.5 kms~^, suggesting an inverse- 
C bisector. The variations of the core intensity up to ±10 % and the core width 
^1 up to ±5 % about the respective averages are much less pronounced, but still 

f--^ . detectable. The core intensity variations lag behind the Dopplershift variations 

fSJ ' about 2.1 min. The Ha core width tends to correlate with the Dopplershift and 

("^ , anticorrelate with the asymmetry, suggesting that more redshifted Ha profiles 

^T ' are wider and the broadening of the Ha core is accompanied with a change 

^^ ' of the core asymmetry from redward to blueward. We also found a striking 

^S , anticorrelation between the core asymmetry and the Dopplershift, suggesting 

fT^ ' a change of the core asymmetry from redward to blueward with an increas- 

ing redshift of the Ha core. The data and the applied analysis do not show 
meaningful tracks of Alfven waves in the selected network element. 
Key words: Sun: chromosphere - Line: profiles 

...' 1. Introduction 

Alfven waves have been invoked as one of possible mechanisms heating the solar 
corona (Alfven, 1947; Osterbrock, 1961). These magnetic transporters of con- 
vective energy are capable of penetrating through the solar atmosphere without 
being reflected and refracted. Since Alfven waves are incompressible, their prop- 
agation through the atmosphere is not associated with density changes seen as 
periodic variations of intensities and line-of-sight velocities (Erdelyi & Fedun, 
2007; Mathioudakis et ai, 2012). Observation of a slanted magnetic fluxtube 
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could reveal Alfven waves as periodic variations of nonthermal broadening of a 
spectral line. This approach was employed in the study by Jess et al. (2009) who 
found periodic oscillations of nonthermal broadening of the Ha spectral line seen 
over a cluster of bright points associated with a distinct upflow without signifi- 
cant periodic variations of intensities and line-of-sight velocities. They interpret 
this as a manifestation of Alfven waves in the solar atmosphere. However, there 
is some ambiguity whether they deal with spectral or integrated intensities. 

Motivated by the results of Jess et al. (2009), we search for signatures of 
Alfven waves in a bright network element occupied by a cluster of G-band 
bright points. We investigate spectral characteristics of the Ha spectral line, 
focusing on temporal variations of the line core width, intensity, Dopplershift, 
and the asymmetry derived by a 4*'^-order-polynomial fitting of a coarse five- 
point sampling of the Ha profile. Since we employ the correcting method and 
results of Koza et al. (2013), and frequently refer to it as Paper I hereafter, we 
invite the reader to have it handy. 



2. Observations 

We use sequences of Ha images of a very quiet area near the disk center (/i = 
0.91) recorded by the Dutch Open Telescope (DOT) on 19 October 2005 from 
09:55:20 UT to 11:05:47UT by the tunable Lyot filter with a FWffMpassband 
of 0.25 A (Gaizauskas, 1976; Bettonvil et al., 2006). The sequence consists of 
scans taken in five wavelengths across the Ha line profile: line center, ±0.35 A 
and ±0.7 A, followed by a single image in the line center. The zero wavelength 
of the filtergrams is centered on the minimum of an average profile taken before 
the observations over the field of view. A time step between the five- wavelength 
scans of the Ha line profile and line center images is 30 s. Since the images at 
two subsequent wavelengths are separated by the 4-to-5-s intervals, one full 5- 
point scan of the Ha profile took about 20-25 s. Within this interval no major 
changes of chromospheric scenery are assumed. The scans were taken in 20- 
frame bursts per particular wavelength setting and then speckle-reconstructed 
applying a Keller - von der Liihe two-channel reconstruction (Keller & von der 
Liihe, 1992). The line center images were obtained in 100-frame bursts restored 
by a full single-channel speckle reconstruction. In this study, we employ 71 Ha 
scans at a regular 60-s cadence taken between 09:55:20 UT and 11:05:22 UT. 
The seeing was only fair, with the Fried parameter ro slightly increasing from 
about 7.5 at the beginning to 9 on average at the end of the observation. As a 
context data, we also use the sequences of 142 G-band and CallH images taken 
synchronously with the Ha images at a regular 30-s cadence. The temporal 
mean of the Ca II H images and examples of the Ha images taken during the 
period of the best seeing are shown in Figs. 1-3. The size of the CallH and Ha 
images is 1112 x 818 px^, with the angular size of pixel of 0.071 arcscc. Details on 
the DOT, its tomographic niultiwavelength imaging, speckle reconstruction, and 
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Figure 1. Images of the quiet-Sun network recorded by DOT on 19 October 2005. Top; 
the temporal mean of the CallH image sequence over its 71-min duration. Bottom: 
the Ha image taken at 10:36:21 UT in the line center at the moment of the best seeing 
occurred in the 41''' min after the start of the observation at 09:55:20 UT. The square 
defines the network element selected for this study and displayed in detail in Figs. 4-7. 
Field of view: 79 x 58 arcsec^ . Tickmark spacing: 2 arcsec. 



J. Koza, P. Siittcrlin, P. Goniory, J. Rybak and A. Kuccra 








n ^'■' ^ 



Figure 2. DOT Ha images at AA = -0.7 A (top) and -0.35 A {bottom) from the line 
center corresponding to the Ha hne-center image in Fig. 1. 
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Figure 3. DOT Hq images at AA = +0.7 A (top) and +0.35 A {bottom) from the line 
center corresponding to the Ha images in Figs. 1 and 2. 
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Figure 4. Temporal evolution of the selected network element seen simultaneously in 
the DOT G-band, CallH, and Hq images taken at five wavelength settings of the Ha 
Lyot filter. The time min corresponds to the moment of the start of the observation 
at 09:55:20 UT. Field of view: 4.3 x 4.3 arcsec . Tickmark spacing: 1 arcsec. 
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Figure 5. The same as in Fig. 4. 
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Figure 6. The same as in Fig. 4. 
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Figure 7. The same as in Fig. 4. 



14 J. Koza, P. Siittcrlin, P. Goniory, J. Rybak and A. Kuccra 



standard reduction procedures are given in Hammerschlag & Bettonvil (1998) 
and Rutten et al. (2004). 

The target area encompasses a full supergranular cell surrounded by a promi- 
nent network as seen in the temporal mean of the CallH sequence in Fig. 1. 
Visual comparison of both panels in Fig. 1 indicates that the brightest network 
elements coincide with roots of long chromospheric fibrils. These dominate Ha 
scene and emanate out from the network elements outlining the structure of 
magnetic field at chromospheric levels. DOT Ha scans render a tomographic 
view of the solar atmosphere from the photosphere (Figs. 2 and 3: ±0.7 A) up 
to the chromosphere sampled by the images in the line center and ±0.35 A off 
(Kontogiannis et a/., 2010). In the Ha +0.7 A image (Fig. 3, top), we can clearly 
see granulation which is missing in the Ha —0.7 A image (Fig. 2, top) due to 
Doppler cancellation. Both images display dark elongated streaks which are 
highly-dopplershifted parts of fibrils seen better in the line center and ±0.35 A 
images. Doppler cancellation enhances the visibility of magnetic elements seen as 
bright points in the Ha —0.7 A image (Fig. 2, top), studied in detail in Leenaarts 
et al. (2006). The square defines the network element selected for this study. The 
size of the square is 61 x 61 px^, corresponding to 4.3 x 4.3arcsec^. 

Figs. 4-7 show enlargements of the square providing a multispectral tomo- 
graphic view on evolution of the selected network element unfolded in time with 
a 4-min timestep from the lower photosphere (G band) up to the chromosphere 
(Ha line center). The G-band cutouts show a cluster of bright points consti- 
tuting the network element. The selected area also involves two large granules, 
with maximum diameters larger than 2 arcsec, occurring in the lower half from 
16 to 28min (Fig. 5) and from 44 to 68min (Figs. 6 and 7). The CallH cutouts 
show brightenings whose overall morphology is often very similar to that of the 
cluster of the G-band bright points. But the former are strongly diffuse because 
of resonance scattering within the solar atmosphere and through expansion of 
idealized magnetostatic fluxtubes with height (Leenaarts et al., 2006). The Ha 
outer wing cutouts (±0.7 A) mix photospheric and chromospheric features be- 
cause of a double-peak contribution function at Ha wings (Schoolman, 1972; 
Leenaarts et al., 2006, 2012). Striking, almost one-to-one correspondence be- 
tween the G-band bright points and the blue wing bright points seen in —0.7 A 
(compare the G-band and Ha —0.7 A cutouts, e.g., at 0, 12, 32, 60, and 64min) 
was explained in Leenaarts et al. (2006). The line center and ±0.35 A cutouts 
bear still signatures of a network brightening, but modulated much with a fast- 
changing blanket of the dark fibrils occurring higher up. 

3. Method, spectral characteristics, and their corrections 

The coarse five-point sampling enables us to construct 1112 x 818 x 71 instan- 
taneous proxy profiles of the Ha spectral line at each pixel. We fit the five- 
wavelength samples of the proxy profiles by a 4*'^-order polynomial to derive 
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the four basic profile measurements: the core intensity Iq, the core velocity vc, 
the core width FW, and the bisector velocity ubi- The values of Iq and vc 
represent the fit minimum and its Dopplcrshift, respectively. The core width is 
the wavelength separation of the two fit flanks at half of the intensity range 
between the fit minimum and the average of the endpoint intensities at ±0.7 A. 
The core width is computed at the average intensity /pw employing the ref- 
erence intensity /rfi defined by equations (3) and (4) in Paper I. Finally, the 
bisector velocity is the velocity equivalent of the wavelength separation between 
the midpoint of the fit at the average intensity /pw and the fit minimum. A 
more detailed definition of the spectral characteristics is given in Paper I. In case 
of a large Dopplcrshift, when the average intensity is larger than the minimum 
of endpoint intensities, i.e., /pw > niin(/_o.7, /+0.7), the fits are extrapolated 
beyond the range ±0.7 A using the polynomial coeflficients and the extrapolated 
fits are used for determining the spectral characteristics. In very rare cases of 
a polynomial undulation, with local maxima lower than the average intensity 
/pWj the algorithm switches to a parabolic fit of the five-point proxy profile to 
determine FW, while Iq and vq are always from the 4"^-order-polynomial fit and 
v-Qi is undefined. The extrapolation (Fig. 8) was invoked just in 0.024% and the 
parabolic fit in 5 x 10""* % out of all pixels. The intensities shown in Figs. 8 and 
9 were normalized with respect to the mean intensities resulting from spatio- 
temporal averaging of the respective image cubes. The five mean intensities were 
converted to the atlas intensity scale convolving the Ha atlas profile with the 
transmission profile of the DOT Ha filter (Paper I, equation (1)). Therefore, the 
integrated intensities in Figs. 8 and 9 are in the wavelength unit. 

The five intensity samples have an instrumental wavelength scale ranging 
from —0.7 to 0.7 A. The wavelengths of particular core fit minima Ac are mea- 
sured with respect to this instrumental wavelength scale with milliangstrom res- 
olution and then converted to the core velocity as: vq = c(Ac — ArJ/Aq , where 
c is the speed of light, Aq is the central wavelength of the Ha line 6562.8 A, and 
the reference wavelength Ar is the spatio-temporal average of all 1112 x 818 x 71 
wavelengths of fit minima Ac over all pixels within the field of view, making no 
distinction between network and internetwork pixels. The reference wavelength 
Ar, is about —18 mA. An adopted sign convention indicates a redshift of the line 
core fit for positive Dopplcrshift and vq . The bisector velocity wri quantifies an 
asymmetry of the line core fit with respect to its minimum. The sign conven- 
tion adopted in Paper I suggests positive vri for a line core fit with a redward 
asymmetry. 

To eliminate the modulation of the spectral characteristics induced by the 
fast-changing blanket of dark fibrils, we construct 71 proxy Ha profiles by spatial 
averaging of intensities within the selected network element (Figs. 4-7) at each 
wavelength. Fig. 9 shows examples of these spatially averaged intensities, their 
polynomial fits, average intensities, and core velocities. Spectral characteristics 
of these spatially averaged intensities are determined solely from their 4"^-order- 
polynomial non-extrapolated fits. To asses an infiucncc of the averaging on the 
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Figure 8. Examples of observed intensities (squares), their 4*'^-order-poIynoiniaI 
(solid), extrapolated (dotted), and parabolic (dashed) fits taken from various loca- 
tions X, y over the selected network element at the indicated time t. The horizontal 
lines indicate the average intensity defining the fit width and the bisector velocity. 
Their solid, dotted, and dashed linestyles correspond to the 4' -order-polynomial, ex- 
trapolated, and parabolic fits, respectively, used in determination of the fit width and 
the bisector velocity. The coordinates x, y refer to the lower left corners of images. 
The selected samples do not represent majority, but rather rare cases with a high core 
velocity vc and highly asymmetric profiles (the middle and bottom panels). 
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Figure 9. Examples of observed intensities (squares) averaged over the selected net- 
work element (Figs. 4-7) and their 4*'^-order-polynomial fits (solid) at the given time 
t. The horizontal lines indicate the average intensity defining the fit width and the 
bisector velocity. The samples illustrate the fits with small (top), medium (middle), 
and high (bottom) redshift with the core velocity vc- 
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final results, we compare them with the instantaneous characteristics obtained 
by the spatial averaging of FW, Iq, vq, and Ubi within the selected area (Figs. 4- 
7). Both, the instantaneous spectral characteristics inferred per pixel and also 
from the spatially averaged intensities are corrected for the deviations estimated 
in Paper I comparing the measurements of FW, Iq, vq, and vbi performed at 
the dopplershifted Ha atlas profile convolved with the filter transmission with 
their reference values. Then the corrected core velocity vq is the solution of the 
equation /(wc) + vq — i> — , where /(wc) is the absolute deviation Av and 
v is the measured core velocity. The corrected FW results from the equation 
FW = fw — AFW{vc) , where fw is the measured fit width and AFW{vc) 
is the absolute deviation AFW, but interpolated on the corrected core velocity 
vc- The corrected ubi is computed in the same way. Due to different intensity 
scales, the corrected Ic results from the equation /c = i/(l + A/r(wc)/100) , 
where i is the measured core intensity and AI^{vc) is the relative deviation A/r, 
but interpolated on the corrected core velocity vc- 



4. Noise, data uncertainties, and references 

Reliability of results inferred by fitting depends critically on the noise presented 
in the data and various intrinsic uncertainties. Therefore, in this section we 
give a brief account of uncertainties of the fitted intensities and the reference 
wavelength defining the core velocity. 

We estimated the upper limit of the noise level Ictnoise presented in the 
Ha filtergrams from la of the temporal variations of the intensity at each pixel 
assuming that Ictnoise < lcr/10. In other words, we assumed that the noise level 
is about one order smaller than la of the temporal variations of the intensity 
due to physical processes on the Sun within the 71-min observing period. This 
rough estimate points out that these Ictnoise errors of the intensities in Fig. 8 
are smaller than the height of the square symbols and can therefore be neglected. 

Similarly, we estimated the noise level of the spatially-averaged intensities 
in Fig. 9 from Ict of the spatial variations of the intensity over the selected area. 
In other words, we assumed that the per-pixel noise level within the area is 
Ict of intensity variations over the area. However, the spatial averaging of the 
intensities reduces their noise level by a factor of 1/^/N or 1/61, where N is the 
number of averaged pixels. This again points out, that the errors of intensities 
in Fig. 9 are much smaller than the height of the square symbols and are thus 
insignificant. 

In this paper, we define the reference wavelength Ar, or zero point of the 
core velocity, as the spatio-temporal average of the wavelengths Ac of the core 
fit minima determined at each pixel. We estimated uncertainty or stability of Ar 
during the 71-min observing period in two ways. First, we constructed 71 proxy 
Ha profiles from instantaneous spatially-averaged intensities, determined the 
wavelengths of their core fit minima, and computed the median of their absolute 
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differences from tlie total average. Second, we took 71 spatial averages of the 
wavelengths of core fit minima for each moment of the observation and computed 
again the median of their absolute differences from the total average. Both 
estimates gave consistent results, indicating that uncertainty of the reference 
wavelength Ar of the core velocity is about ±0.3 km s^^. 

When determining the reference wavelength Ar we made no distinction be- 
tween the network and internetwork areas and this might be a source of a 
systematic offset. We checked this feature by selecting only large internetwork 
areas for determination of Ar. This one is blueshifted about 0.2kms~^ with 
respect to the reference wavelength from the whole field of view. 

5. Results 

This section presents a pictorial overview of results inferred from the data un- 
dergone various processing. Gray shades in Figs. 10 - 12 highlight differences 
between the uncorrected (gray) and corrected (black) data, while the differ- 
ent styles distinguish the results of intensity averaging (solid lines) and spatial 
averaging of spectral characteristics (symbols) over the network element. 

5.1. Temporal evolution of the spectral characteristics 

The top panel of Fig. 10 shows temporal evolutions of the relative variations 
of the fit width FW and the core intensity Iq of the line core fits of the Ha 
spectral line in the network element (Figs. 4-7) with respect to their temporal 
means. The bottom panel shows temporal evolutions of the core velocity vq and 
the bisector velocity vbi, whose positive values indicate redshift and redward 
asymmetry of the line core. All these characteristics were derived from the orig- 
inal Ha datacubes downloaded from the DOT database^ . An application of the 
corrections (Paper I) increases apparently the amplitude of the FW^ variations 
after the 40"^ min, but decreases both velocities within the whole 71-min pe- 
riod. The values of FW and 7c seem to be largely insensitive to the choice of 
averaging for the size of the selected area, but this is not the case of vc and wbi. 
Their corrected values resulting from spatial averaging over the network element 
(black diamonds and triangles in the bottom panel) differ systematically from 
the values inferred from spatial means of the intensities (black thick and thin 
lines). This is probably due to high non-linearity of the respective correcting 
functions shown in Paper 1. The temporal evolutions of the characteristics do 
not show any apparent systematic trend and fast short-period variations dom- 
inate them. The variations of corrected FW are within the range ±5% around 
the mean. The variations of Iq indicate a rapid darkening of the selected net- 
work element (Fig. 4) from 11 % above an average to average within the first 10 
minutes. Later on, Iq variations lie within the range from —9% to +6%. The 
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Figure 10. Temporal evolution of the fit width FW, the core intensity 7c, the core 
velocity vc , and the bisector velocity wbi of the line core fit of the Ha spectral line in 
the network element (Figs. 4-7) derived by the 4"'-order polynomial. FW and /c are 
shown as variations with respect to their temporal means (top panel) . Positive vc and 
vbi signalize a redshift (middle and bottom panel). Gray and black lines indicate the 
uncorrected and corrected characteristics, respectively, derived from spatial means of 
the intensities. Gray and black symbols shown with the time step of 2 min indicate 
the uncorrected and corrected characteristics, respectively, derived as spatial means of 
particular values at each pixel. 



corrected vc displays redshifts up to 4knis~^, with an average of 1.5 kms""'^ 
and a noticeable correlation with the corrected FW, most apparent after the 
fortieth minute. In the same period, a large expanding granule with the maxi- 
mum diameter larger than 2 arcsec appeared in the low photosphere (see Figs. 6 
and 7). Fig. 10 suggests also some similarity between the vq and Iq variations. 
The Fourier cross correlation indicates that Iq lags behind vq about 2.1 min 
on average. The Ha core possesses mostly a redward asymmetry corresponding 
to an inverse-C bisector with the average bisector velocity of about 0.5 kms^^ 
with excursions up to Iknis^^. 



Search for Alfvcn waves in a bright network clement observed in Hct 



21 





1.06 






1.04 


<> w <§> 


xC 


1.02 


_ ,^ ':-^<^ ^ o o^ 


^ 


1.00 


o ^ o 




0.98 


o o 




0.% 


. . ^1 1 





1.06 


1 1 1 { 1 1 1 




1.04 




<. 


1.02 




t 


1.00 


- ^?0 CC = 0.24 _ 
^ ^ CC = 0.71 




0.98 








O.W) 


1 O ,. 1 , t 



5200 



6400 - 



5600 fiWX) 



2 4 

Vf. [km s" ] 



6000 



5600 



5200 






' 1 


1 'O 1 ' 




(i40(J 


— 


-') 


^ 




^ 




■ 




o 

- CO ^ 


o^to 


- 


6ffl» 


^ 


- 




'- <%^ 


^^ o^ 


; 




<>^ 


^ "^ o 


■ 


5600 


^ 


^ o 


- 






^ 


■ 


5700 


1 


1,1, 





2 4 

V(^ [km s' [ 



2 ~r 



^^^^ 




o 



ofM 



^#v 



cc = .().! 7 O 
_ CC = -0.42 O O _ 
J I I I I L 



1 - 



..,. . I. 


' 1 ' [ ' 
o 


■ 


y81 


^^o^ 


- 


- 


'%. 


- 


. cc = 


= -0.93 1>3 ^ 


■ 


_ cc = 

1 


= -0.68 OO 
1,1, 


- 



0.96 0.98 I.OO 1.02 1.04 1.06 
FW [A[ 



2 4 



Figure 11. Scatter plots of the spectral characteristics shown in Fig. 10 - the case of 
spatial means of the intensities. Gray and black symbols indicate the uncorrected and 
corrected characteristics, respectively. Unlike Fig. 10, the plots show absolute values of 
F IF and Ic- Positive vc and «bi signalize a redshift. There are also shown the corre- 
lation coefficients CC for the pairs FW - vc, i^bi - FW, and «bi - ^"0 for uncorrected 
(gray) and the corrected (black) values. 
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Figure 12. Scatter plots of the spectral characteristics shown in Fig. 10 - the case 
of spatial means of the spectral characteristics at each pixel. Gray and black symbols 
indicate the uncorrected and corrected characteristics, respectively. Unlike Fig. 10, the 
plots show absolute values of FW and Ic- Positive vc and ubi signalize a redshift. 
There are also shown the correlation coefficients CC for the pairs FW- vc, wbi - FW, 
and uei - ^c for the uncorrected (gray) and corrected (black) values. 
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5.2. Scatter plots 

The spectral characteristics from Fig. 10 are displayed in the form of scatter 
plots in Figs. 11 and 12, showing separately the results of spatial averaging of 
the intensities and the spectral characteristics at each pixel, respectively. The 
top left panels show FW against /c which do not correlate, neither uncorrected 
nor corrected. The top right panels suggest a positive correlation between FW 
and vc, which occurs only for the corrected values. The correlation implies that 
more redshifted profiles tend to be wider. The middle panels indicate that there 
is no apparent correspondence between Iq and the velocities. The values of Iq 
do not correlate with vq, not even after shifting Iq forward (i.e., to the left with 
respect to the zero of the time axis) about 2.1 min to compensate its apparent lag 
indicated in Fig. 10. The bottom left panel of Fig. 12 suggests an anticorrelation 
between wbi and FW, which occurs only for the corrected values resulting from 
their spatial averaging at each pixel. The anticorrelation implies that while 
narrower Ha cores display a redward asymmetry, the wider and more redshifted 
profiles (see the FW- vc correlation above) display a blueward asymmetry. The 
bottom right panels also suggest a striking anticorrelation between wbi and vc ■ 
Notably, the profiles with the zero Dopplershift display a redward asymmetry 
with wbi = Ikms""'^ and become symmetric for vc of 2kms~^, but tend to 
display an increasing blueward asymmetry with further-increasing redshift. 

6. Discussion 

We derived spectral characteristics of the Ha spectral line observed in the net- 
work element with the aim of seeking signatures of the Alfven waves as reported 
in Jess et at (2009). How do our findings compare to the results of these authors? 
While Jess et al. (2009) reported an average blueshift of 23kms~^ in a cluster of 
bright points seen in Ha, we found an average redshift of about 1.5 km s^^ . Obvi- 
ously, this is a large discrepancy and we suspect that our particular choice of the 
reference wavelength of the core velocity might be partially responsible for that. 
We remind that we made no distinction between the network and internetwork 
(see Sections). But as we showed in Section 4, selecting only internetwork areas 
shifts the reference wavelength blueward only insignificantly up to 0.2kms~^ 
and the value of the average redshift is significantly above the uncertainty of 
0.3kms~^ of the reference wavelength. Further, Jess et al. (2009) reported an 
absence of any significant intensity variations in the cluster of bright points. 
However, there is some uncertainty whether they refer to the spectral intensity 
in the line core, or the integrated intensity. Our results suggest that the core 
width variations in the selected network element are associated with both the 
core velocity and core intensity variations. Finally, the authors detected FWHM 
oscillations of the Ha spectral line with the strongest power in the 400-to-500 s 
interval. We postpone detailed cross correlation, Fourier, and wavelet analyzes 
of our results to a forthcoming paper. We conclude that, most likely, a differ- 
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ent mechanism works in the selected network element than in the cluster of 
the bright points studied by Jess et al. (2009). The inversc-C bisector and the 
average redshift of the Ha line core are more symptomatic for propagation of 
chromospheric shocks, as suggested in Uitenbroek (2006) and demonstrated in 
numerical simulations by Heggland et al. (2011). Then the lagging of the core 
intensity /c behind the core velocity vq variations about 2.1 min might be char- 
acteristic either for the p-mode oscillations leaking into the chromosphere, or for 
the shock wave propagation. Although there is an apparent similarity of their 
variations best seen after the fortieth minute (Fig. 10), shifting /c forward about 
the 2-min lag does not improve the Iq - vq correlations shown in the middle 
right panels of the scatter plots in Figs. 11 and 12, which may be ascribed to a 
variable phase shift which varies in time and/or frequency. 

Since the standard DOT observations do not involve any precise wavelength 
calibration of the Ha filtergrams, we refrained from interpreting straightfor- 
wardly redshifts and blueshifts as actual mass downflows and upflows, respec- 
tively. Despite these facts, we will try to reason that most of the positive values 
of vq indicated by the redshifts of the Ha core may represent the actual mass 
downflows. As we noted in Section 2, the zero wavelength of the filtergrams is 
centered on the minimum of an average profile taken before the observations 
over the field of view. Thus, the zero wavelength approximately corresponds 
to the Ha solar disc-center wavelength of Asoi = 6562.81 A according to the 
spectral FTS atlas by Neckel (1999), whose wavelength scale was already cor- 
rected for the rotational and radial velocity of the Earth. Consequently, the 
reference wavelength given in Section 3 corresponds to Ar = Agoi — 0.018 A — 
6562.79 A. The local standard of rest at the solar surface for the Ha spectral 
line is: Aq — Aiab(l + ^gr/c) = 6562.80 A, where Aiab = 6562.79 A is the air 
laboratory wavelength of Ha adopted from the NIST database (Kramida et. al, 
2012), WGR = 0.636 km s^^ is the solar gravitational redshift, and c is the speed 
of light. Then the core velocities of Ha measured within this study are offset 
only about (Aq — Ar)c/Ao = 0.5 km s^^ with respect to the local standard of 
rest at the target area, with the following consequences: 

— the zero of vq and all values of uri in Fig. 10 should be shifted redward (i.e., 
up) about 0.5kms^^, 

— the zeros of the vc axes in Figs. 11 and 12 should be shifted redward (i.e., 
to the right) about 0.5kms~^, 

— most of the positive values of vq indicated by the redshifts of the Ha core 
may represent the actual mass downflows. 

7. Conclusions 

A selected network element exhibits distinct spectral characteristics seen in their 
temporal evolutions and scatter plots correlating pairs of values from the same 
instant. In general, this involves: 
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— the fit width variations about ±5%, i.e., a few tens of milliangstroms with 
respect to the average of 1.01 A, 

— the core intensity variations about ±10% with respect to the average, 

— the core velocity variations from to 4knis^^ about the averageof 1.5kms~^, 
implying a redshift, 

— the core asymmetry variations of the Ha hne core up to 1 knis^^ about the 
average of 0.5kms~^, suggesting an inverse-C bisector. 

The Ha core width tends to correlate with the Dopplershift and anticorrelate 
with the asymmetry, suggesting that more rcdshifted Ha profiles are wider and 
the broadening of the Ha core is accompanied with a change of the core asymme- 
try from redward to blueward. We found also a striking anticorrelation between 
the core asymmetry and the Dopplershift, suggesting also a change of the core 
asymmetry from redward to blueward, with an increasing redshift of the Ha 
core. A question is which of these patterns are unique just for the network and 
which are also common in the internetwork. This can be resolved in a stand-alone 
study, comparing the spectral characteristics of the Ha line in the network and 
internetwork. Finally, an absence of blucshift and detected intensity and velocity 
variations probably exclude a presence of Alfvcn waves in the selected network 
element according to the criteria given in Jess et at (2009). 
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